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ABSTRACT 

For high precise baseline and source position determinations with VLBI, large bandwidth techniques 
are required to achieve the necessary delay resolution of < 0.1 ns. In local interferometry, the 
interferometric phase observable offers an equally high resolution in spite of the much smaller 
bandwidths used. In the present paper, the conditions for the applicability of the phase observable 
in very long baseline interferometry (VLBI) are discussed and examples of adequate observation 
programs are given. These considerations apply in particular to the Mark II system, which is already 
in operation at many observatories. 
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DELAY AND PHASE OBSERVATIONS 


In radio interferometry the travelling time r of a wavefront between the antennas at two sites is 
measured by maximizing the cross-correlation function of the two signal streams received at the two 
antennas. Simultaneously, the phase and the rate of the interference fringes are determined in the 
processor as they vary with the Earth’s rotation. During correlation, the data stream from one sta- 
tion is delayed quasi-continuously in such a way that the changing geometric delay r s is almost 
completely compensated for. This gives rise to a rather low residual fringe frequency the phase of 
which slowly varies on the scale of a few turns per minute. In order to be able to observe those 
fringes which have maximum amplitude, the correlation is carried out simultaneously in a number 
of delay channels separated by 


At 


1 

2 B 0 
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where B 0 is the effective bandwidth of the processed signal streams. The output of the correlator is 
usually described by the complex cross correlation functions, which translates the response of an 
interferometer system to a point source*: 


R (t, r A ) = W SNR • D (r - tj) ■ e 
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where 

W SNR 

t 


a factor depending essentially on the signal-to-noise ratio (SNR) of the interferom- 
eter system observing a given source, 
observing time (i.e., UTC), 

actual delay including all instrumental and atmospheric effects, 
discrete time delays as computed for the different delay channels. 


D (r - Ti ) 
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7T B q (r - Tj) 


D = delay resolution function, 


(2a) 


0f (t, Tj) = residual fringe phase (including all instrumental and atmospheric effects), 


and 


<j)f = * t 


(2b) 


where cof = residual fringe frequency. 


(For more details the reader is referred to Thomas, 1972 or Moran, 1976.) 


For the discussion that follows, it is important to note that the interferometric observables used in 
geodetic VLBI are all contained in equation (2), namely the group delay r(t), the phase 0f (t), and 
the fringe frequency cof (t). 


♦The source structure is expressed by an extra amplitude and phase term which is omitted here for simplicity. 
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The group delay is estimated from the delay resolution function (2a), which becomes a maximum 
for Tj = r. The accuracy with which the exact fraction of Tj can be determined, depends mainly on 
the halfwidth of the delay function, which is given by 

Ar h =4"' (3) 

Bo 

There are different possible methods of delay estimation* (i.e., Whitney et al., 1976) all of which 
yield a precision a T of roughly 1 percent of Ar^, depending on the SNR and the available integra- 
tion time per observation. The following table is meant to give an overall idea on the delay estima- 
tion performance of different VLBI systems. 


Table 1 


System 

B Q [MHz] 

Ar h [ns] 

a T [ns] 

c • a T [cm] 

Canadian analog. 

4 

250 

2.5 

75 

Mk I 1 channel 

0.36 

2,778 

28 

840 

Mk I synth. 8 channels 

23** 

43 

0.5 

15 

Mk II 1 channel 

2 

500 

5 

150 

Mk III 28 channels 

56** 

18 

0.2 

6 

Sat-Link (concept.) 

100 

10 

0.1 

3 


* *high-degree of flexibility (Mk III: 0.2 - 300 MHz!) 


The instrumental and atmospheric errors are discussed in connection with the phase observable. 

The residual fringe phase and the fringe rate are obtained from the sine and cosine parts of the 
cross-correlation function (last term in equation (2)). Due to the close relationship of af and cuf 
(equation (2b)), these observables are determined simultaneously, either from an ordinary sine wave 
adjustment or using the fourier transform into the frequency domain, where the maximum of S 
(cufj) is estimated. After convolution with e'J^f*, the phase for a given instant t can be computed 
from the real and imaginary components at lag zero. These methods, which allow to establish the 
function 0 (t) over a certain interval of time (usually the duration of an uninterrupted source scan), 
are often referred to as “phase tracking” methods (Thomas, 1972). The accuracy of phase tracking, 
which is inherent to the system, only weakly depends on system bandwidth through SNR consider- 
ations. So in the case of observations on strong sources, no marked differences in accuracy can be 
found comparing the different VLBI systems. The tracked phase can be converted into a so-called 
phase delay by using 

tM) = — — 0 (t) + 0 O (4) 

v 2 * v o 


Rigorous delay estimation is done iteratively together with the fringe rate estimation: D = D(r, OJf). 
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where v 0 = center of observing frequency band, and 

(j> 0 - constant phase zero term, accounting for the unknown number of turns at the start 
of each source scan. 

The phase delay accuracy, which is equivalent to the phase accuracy, is composed of the small phase 
tracking error plus the much larger instrumental and atmospheric phase errors listed below: 

— clock stability (phase variations of the local oscillators) 

— instrumental delay (electronics, cables + mechanical structure) 

— atmospheric delay (ionosphere; troposphere, dry and wet) 

— resolved sources add phase variations due to their structure. 

Typical phase tracking errors on strong sources are about one-tenth of a turn or smaller; this trans- 
lates into phase delay errors of some tens of picoseconds only, corresponding to the observing fre- 
quency used. The overall phase variations due to the above mentioned effects, however, amount to 
about 0.1 ns (if hydrogen masers are used as frequency standards). It is difficult to give a represen- 
tative number here because of the wide variability of atmospheric and instrumental conditions. It 
may be possible in the near future to reduce this number by a factor of two or three using satellite 
LO links, instrumental delay calibration, dual-frequency receiving systems and water vapor radiom- 
etry. But even in the less favourable case of a T ~ 0.1 ns, the phase delay accuracy of any narrow 
band system is equivalent to a group delay precision of a wide bandwidth system with B 0 = 100 
MHz. 

What can be done to take advantage of the phase observables for meaningful geodetic work, will be 
shown in the next paragraphs. 


GEODETIC BASELINE DETERMINATIONS USING PHASE OBSERVATIONS 

There are several ways in which phase observations are made following the purpose of the experi- 
ment. 


Method A: Uninterrupted Scans for Mapping the Structure of Compact Radio Sources 

In most of these experiments the telescopes are aimed at the same source for many hours, with the 
exception of a few short calibration scans. During each of these scans, the phase can be tracked 
without ambiguity problems, even at higher frequencies, as long as interference fringes are seen. 

For a unique baseline solution with this kind of data, the following conditions have to be met: 

1 . The polar component (z = component) of the baseline has to be kept fixed. 

2. At the start of each source scan, a new phase zero unknown (which is equal to a clock offset 
unknown) has to be introduced. 
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3. The duration of the sources scans should not be shorter than 3 to 4 hours. An optimum is 
reached if the hour angle range covers 180° (=12 hours). 

These conditions can be deduced from the observation equations, where the baseline terms ex- 
pressed in cylindrical coordinates are 


dbi - 


1 


1 - — sin 6 

1 c 

_ 1 


(polar component) 


db 2 = — cos 6 cos h (equatorial component) 


dAu = _Acos5sinh (longitude of baseline vector) 
c 

with 5 = declination and h = hour angle of the observed source. 


An example of a baseline solution with 13 cm data on an intercontinnetal baseline is shown in fig- 
ure 1 . These data are the first that could be used with a program completed only very recently. In 
spite of the poor quality of these data (rubidium clocks were used at both stations), a comparison 
with the delay solution of the same experiment demonstrates the substantial gain in accuracy: 


delay solution: 

dbj 

= 47.52 

± 

1.62 m 

o T = ±9.5 ns 

db 2 

= 16.00 

+ 

2.86 m 


d\ b 

= 4797 

+ 

0732 

phase delay solution: 

dbj 

= 



a T = ±0.20 ns 

1 <P 

db 2 

= 19.32 

+ 

0.24 m 

d\ b 

= 4742 

+ 

0703 


Another example is shown in figure 2, where phase data on the 832 km baseline Effelsberg-Onsala 
at 1.3 cm (22.2 GHz) were converted into phase delays. Due to the small ambiguity spacing at this 
high frequency, it proved to be difficult to connect the phases across data gaps larger than 5 minutes 
(coherence time of the interferometer). 

The phase delay solution with scan on 3C273 and 3C84 gave the following results: 

a T = ±0.15 ns db 2 = 0.15 ± 0.054 m 

d\ b = 3736 ± 07008 

The drawback of method A, however, is twofold: the baseline solution is incomplete and only a 
few sources can be observed during one experiment. Moreover, very severe bounds have to be 
placed on the systematic effects due to the instrumentation and the atmosphere because there is a 
high sensitivity of these solutions to any non-linear effects in the phase delay observables. 
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Figure 1. Phase delay residuals at 2.3 GHz (13 cm! on the baseline Effelsberg-Johannesburg (8024 km). Rb-stand. at Effb 



recidupil: oso-^p --lffelg 
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Figure 2. Phase delay residuals at 22.2 GHz {1 .3 cm) on the baseline Effelsberg-Onsala {832 km). H-masers 
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Method B: Interrupted Scans on Source Pairs (Source Switching Between Pairs) for Accurate Deter- 
mination of Relative Source Positions* 

In this case, conditions for a unique baseline solution are: 

1 . Polar baseline component is fixed. 

2. A phase zero unknown for each observed source has to be introduced. 

3. Same duration minimum as method A, the optimum being again 12 hours. 

4. The length of the individual scans (including telescope drive time) should not exceed a limit set 
by the possibility to interpolate the phase change up to the correct number of turns between 
two consecutive scans on the same source. This condition is satisfied if the integrated phase 
error is substantially smaller than 180°. Remaining phase ambiguities within the scans pertain- 
ing to one source can be removed in an iterative process. A portion of a plot still containing 
some ambiguity errors is shown in figure 3. 

The number of sources that can be observed in this scheme is larger by a factor of two, which in- 
creases the stability of the solution. Some comparative computations based on a constant observa- 
tional error of o T =0.1 ns are listed in the following table: 

0 


Table 2 

Baseline Effelsberg — Green Bank (6 335 km) 


Source 

Pair 

Deck 
[ °] 

Total Observation 
Time 

Baseline Component 
Accuracy 

b 2 


VRO 
3C 371 

42° 

70° 

~7 h 

± 8 cm 

±07015 

CTA 102 
3C 371 

11° 

70° 

~7 h 

± 10 cm 

±070.17 

VRO 
CTA 102 

42° 

11° 

~7 h 

±13 cm 

±07004 


(duration of individual source scans: 24 min) 


♦Shapiro, Wittels, et al., 1979. 
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Figure 3. Phase delay residuals at 2.3 GHz (Effb.-Joh.). Source Switching Mode (3C273 and 3C279) 
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Method C: Interrupted Scans on Many Different Sources 

The sampling of sources all over the commonly visible portion of the sky would, of course, provide 
the most stable least squares solution, a fact which can best be demonstrated by comparing the cor- 
relation matrices of methods A, B and C: unless the observation period is extended over 12 hours 
or more, the solutions of methods A and B tend to show strong correlations (0.9) between some of 
the estimated parameters. It is interesting, however, that the equatorial baseline projection b 2 in all 
cases correlates only very weakly with the other parameters. 

The realization of method C becomes increasingly difficult with an increasing number of sources 
involved in the switching scheme. This is easily understood because the gaps between scans on the 
same source widen considerably, preventing correct phase interpolation. In this case, more sophis- 
ticated methods of ambiguity elimination are required. 


AMBIGUITY ELIMINATION 

In order to obtain a complete baseline solution including the z-component, the remaining phase 
ambiguities between scans to different sources must be resolved. In Rogers et al., 1978 an interest- 
ing method of phase ambiguity elimination is described, which relies on a preliminary wide band 
delay solution. At the relatively high observing frequency of 7.85 GHz, one phase turn corresponds 
to an ambiguity spacing of 3.8 cm or 0.13 ns. This was roughly equivalent to the uncertainty of the 
adjusted group delay observations from the preliminary solution. Now these adjusted delay obser- 
vations were converted to reference phases which were then subtracted from the observed phases. 
After suppression of clock drift effects, the phase differences were treated in a special function 

f(b) = 2 cos A (1)* 

i 

If f (?) is maximized by varying the baseline solution parameters, all the different A0j tend to crop 
around zero phase. The values obtained for b in this way are used to remove the ambiguities from 
the original phase observations, which are then converted to phase delays and processed in the final 
solution. The accuracy gain proved to be considerable: the postfit rms of the phase delay residuals 
was down to about 10 ps compared to some 200 ps** with the group delay solution. 

From the above experiences, which were obtained with the Mark I system on a very short baseline 
(1.24 km), we may draw one important conclusion (which is independent of baseline length): it is 
possible to resolve phase ambiguities if a preliminary solution with a delay phase accuracy of about 
one turn (or even somewhat more) can be achieved. It is obvious that the requirements for the pre- 
liminary delay solution will become less stringent with decreasing frequency. However, a lower 


*This sum was weighted if the observations are not equally spaced. 

**The synthesized bandwidth was ~100 MHz. 
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limit of the observing frequency is set by the ionospheric path length variations, which increase with 
l/i> 0 2 . Experiments are planned by our group to test the available frequencies at the lower end of 
the GHz-scale. 

With reference to the limited Mark II system’s performance, the following configuration appears to 
hold promise without showing much complexity: A basic observing frequency could be chosen at 
1.66 GHz (18 cm) with an ambiguity spacing of 0.6 ns. In order to resolve this ambiguity, a band- 
width of about 10 MHz would have to be synthesized. This could be achieved with four channels 
spread over the band (see figure 4). In order to reduce SNR losses, two recording sets are used in 
parallel, each switching rapidly between two channels. The resulting delay ambiguity of about 300 
ns can be solved with the coarse delay of each of the channels, while the fine delay itself will have a 
halfwidth of ±100 ns. The accuracy of the preliminary delay solution of this configuration should 
be about 0.8 ns with good SNR, which is in agreement with the above mentioned requirements. 

The phase delay accuracy of the proposed system will be limited by the same effects as considered 
in paragraph 1 , this amounts to a few cm or about 0.1 ns on intercontinental baselines. 

The proposed Mark II phase delay system could be equivalent to an economy version of the Mark III 
system with the advantage of being at a lower cost level. On the other hand, if used in conjunction 
with the highly flexible Mark III system, the proposed scheme should result in a substantial reduc- 
tion of tape consumption because only two channels would have to be recorded simultaneously 
(factor of 14). 

The 38 MHz Mark II system (bandwidth synthesis), as it is used at present by the JPL-Caltech group, 
should be particularly suitable for phase delay analysis with ambiguity elimination. 


CONCLUSIONS 

There is no doubt that with wide band group delay techniques the optimal geodetic baseline results 
can be achieved, but it is also clear that these techniques, due to their high costs and complexity, 
will not be in reach of many observatories. 

The test computations with real and simulated data presented here are meant to demonstrate that 
meaningful geodetic work can be done with the existing Mark II technology if the phase observable 
is used. In spite of the incomplete baseline solution (i.e., the z-component cannot be determined), 
it is possible, for example, to measure Earth tides and continental drift on east-west baselines if the 
changes in the equatorial baseline component are considered. Furthermore, the longitude of the 
baseline vector can be used for the determination of UT 1 variations. 

An additional advantage of the proposed phase observation scheme resides in the fact that source 
structure analysis is possible using the closure phase method on the same data. Thus, the structure 
is monitored simultaneously and yet independently of the baseline solution. However, due to the 
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center at 1.66* GHz 



Mk U rec. H 


Mk ”5 rec. B 


a) 10 MHz synthesis 



b) fine delay function 


Figure 4. System characteristics for ambiguity elimination at 18 cm . 
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sensitivity of the pure phase solutions to instrumental and atmospheric effects, the full benefit of 
the high resolution of the phase observable can only be obtained if these effects are monitored with 
great care. 

As an alternative, a system based on Mark II technology is outlined, which is capable of (a posteriori) 
resolving the phase ambiguities in order to achieve a complete baseline solution. With this system, 
it is possible to optimize observing strategies to obtain well conditioned baseline solutions that are 
less sensitive to long term instrumental and atmospheric effects. 


REFERENCES 

Moran, J. M.: Very Long Baseline Interferometric Observations and Data Reduction. Methods of 
Experimental Physics, Vol. 12, Part C, Radio Observations, Academic Press, New York, 1976. 

Rogers, A. E. E. et al. : Geodesy by Radio Interferometry: Determination of a 1.24-km Baseline 
Vector with ~5 mm Repeatability. J. Geophys. Res., Vol. 83, p. 325, 1978. 

Shapiro, I. I., Wittels, J. J., et al.: Submilliarcsecond Astrometry via VLBI: I Relative Position of 
the Radio Sources 3C345 and NRAO 512. Preprint submitted to The Astronomical Journal 
April 1979. 

Thomas, J. B.: An Analysis of Long Baseline Radio Interferometry. The Deep Space Network Pro- 
gress Report, Technical Report 32-1526, Vol. VII, VIII, XVI, Jet Propulsion Laboratory, 
Pasadena, Calif., 1972. 

Whitney, A. R. et al.: A Very-Long-Baseline Interferometer System for Geodetic Applications. 
Radio Science, Vol. 1 1 , p. 421-432, 1976. 


